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By a refinement of the procedure used in the usual earth-flattening approximation, the 
problem of propagation around a spherical earth is reduced to an exact equation of the same 
form. Thereby the earth-flattening procedure becomes applicable to arbitrarily large 
heights and distances. It is also found that existing solutions of the approximate equations 
can be re-evaluated to yield the exact solutions for slightly different refractive index distribu- 
tions. 



1. Introduction 

The problem of propagation of radio waves around 
a spherical earth has been attacked by many writers 
over a period of more than a half a century. For the 
idealized case of homogeneous atmosphere and earth, 
the presently accepted solution is due to van der Pol 
and Bremmer [1]. This makes use of a modified form 
of the Watson transformation to convert from a very 
slowly converging series of spherical harmonics to a 
more rapidly converging series of residue waves. 
Physically, the transformation is from waves pro- 
gressing radially out from the earth to waves which 
progress along the earth's surface. 

In order to reduce the mathematical complexities 
accompanying the van der Pol — Bremmer treatment, 
M. H. L. Pryce [2, 3] introduced the so-called earth- 
flattening approximation, in which plausible approxi- 
mations are introduced to simplify the original dif- 
ferential equations. This procedure reduces the 
problem of propagation over spherical earth to 
propagation over a plane earth with an atmosphere 
having a modified refractive index. In this way the 
van der Pol — Bremmer type of result for the homoge- 
neous atmosphere can be obtained in a much more 
direct manner. Furthermore, the method can be 
applied to inhomogeneous atmospheres in which the 
refractive index varies only with height (so-called 
horizontal stratification) . 

Pekeris [4] investigated the accuracy of the earth- 
flattening approximation in the case of the standard 
atmosphere (in which the modified refractive index 
increases linearly with height). He found that the 
accuracy decreased with increasing distance along 
the earth's surface and also with increasing height 
above the surface. If r is the great-circle distance 
and a the earth's radius, the error with distance is 
proportional to (r/a) 2 independently of the frequency, 
amounting to about 8 percent at r/a=l. The error 
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with height, however, is proportional to frequency, 
and to height as h 5/2 . The height error becomes 
serious in the microwave region at heights in the 
order of a few thousand feet. 

In this paper we improve and extend the Pryce and 
Pekeris treatments to obtain exact differential equa- 
tions for the spherical geometry in terms of equations 
of plane earth type. The solutions of these hold for 
arbitrary heights and distances. Comparison with 
the previous earth-flattening approximation leads to 
an interpretation of the errors therein. 

2. Formulation of the Problem 

We consider propagation from a source above a 
spherical earth of radius a with an atmosphere whose 
refractive index n is a function only of the radial 
distance (i.e., horizontally-stratified atmosphere). 
The treatment is simplified without any essential 
restriction by assuming that radiation from the source 
is independent of azimuth (we shall indicate later the 
extension to remove this restriction). Then it can 
be shown [5] that the electric and magnetic fields at 
any point may be derived from a radial Hertz vector 
\f/ } which satisfies the reduced wave equation 



vV+W=o, 



(2.1) 



where ko=27rl\=2irf/c. By adopting a spherical co- 
ordinate system (i?,0,<£) with origin at the center of 
the earth and with the source on the polar axis, (2.1) 
(with the assumed azimuth independence of \p) 
reduces to 



irdi, 



^)+wk-^^>^-°- 



(2.2) 



This equation can be separated in the usual way by 
letting 

*=T(e)U(R), 
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whereby there result the two ordinary differential 
equations 



d 2 T , 



JT 



de^ +cote de +(as)2T=0 > 



(2.4) 



in which (as) 2 is the separation constant and S=s/k . 
We now introduce the transformations 3 

1=alog(R/a), (2.5) 

r=ad. (2.6) 

Equations (2.3) and (2.4) become, respectively, 

W + l ^+^(^/«-SW=0, (2.7) 



d 2 T 1 dT 

^4+i cot (r/a) ^+s 2 T=Q. 

dr 2 a ' J dr 



(2.8) 



S 2 is the eigenvalue of the differential eq (2.7), which 
eventually is determined by the boundary conditions. 
The problem of solving (2.1) thus is equivalent to 
solving (2.7) and (2.8). 

3. The Height-Gain Function U 

We consider first (2.7) for the height-gain function 
U. By introducing a new variable V, given by 



jj=(a/RyV=Ve-i /{2a \ 



(2.7) becomes 



dW 
dr) 



2+ kl(N*-a°)V=0, 



(3.1) 



(3.2) 



in which we introduce as our definition of the modified 
refractive index 



and define 



N=ne 1,/a =nR/a, 



a 2 =S 2 +l/(2k a) 2 



(3.2) is exactly the same inform as the plane-earth 
equation, which is 



^-°+^(7V 2 -S 2 )f7o=0, 



(3.3) 



where h=R -a is the actual height above the earth's 
surface, and 

N 2 =n 2 (l + h/ay 2 ~n 2 (l + 2h/a). 



3 These are the transformations used by Pryce [3], who credited Pekeris for the 
range transformation and E. T. Copson for the height transformation. 



The method of the usual earth-flattening approxi- 
mation is to approximate U by U , which is the 
solution of (3.3). But U can be obtained exactly 
by means of (3.1), in which V is the solution of (3.2). 
This fact will permit us to interpret properly the 
significance of the earth-flattening approximation 
with respect to the height dependence. 

4. Interpretation of the Earth-Flattening 
Approximation 

Equations (3.2) and (3.3) are identical only in 
form. However, the height variable rj differs from 
h, and, for a given atmosphere, N as a function of 
rj will not have the same analytical form as N as a 
function of h. In addition, the eigenvalue a 2 differs 
from S 2 by the small real quantity l/(2k a) 2 . Thus 
it follows that a given solution U of the approximate 
earth-flattened eq (3.3) corresponds (but is not 
equal) to the exact solution of a slightly different 
problem. The approximate solution differs from 
the corresponding exact solution in three respects, 
namely: (a) By the factor ^ (a/R) 1/2 in (3.1) (this 
factor usually is nearly unity) ; (b) the modified 
index profile is different because of the difference in 
height scale; 4 (c) the eigenvalue is shifted slightly. 
Both (b) and (c) lead to differences between Uo and 
U which involve the frequency, (b) probably is the 
more important effect, so that the error of the earth- 
flattening approximation as a function of height prob- 
ably is due primarily to the choice of height scale. 

5. The Standard Atmosphere 

The standard atmosphere is now defined by 

N 2 (h)=l+gh ) 

ff=3/(2e0. 

In this case, (3.3) becomes Stokes' equation, whose 
solutions are Airy functions, or the modified Hankel 
functions of order one-third [6]. In view of the fact 
that an exact solution is available in terms of our 
function U, it appears that it may be desirable to 
define the standard atmosphere by the relation 

2^(10=1+01. 

In this case (3.2) likewise becomes Stokes' equation. 

6. The Distance Function T 

The distance depend (Mice is given by the function 
T, which satisfies the Legendre differential eq (2.8). 
In the earth-flattening approximation, T is approxi- 
mated by a function T which satisfies the differential 
equation 

d 2 T , 1 dT 



dr 2 r dr 



'+s'T a =0, 



(6.1) 



4 This shall be discussed in greater detail for specific cases in a later paper. 
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which is the Bessel equation of order zero. Equa- 
tion (6.1) gives the distance dependence for plane 
earth. 

In investigating the accuracy of this approxima- 
tion, Pekeris [4] assumed that T could be written 
in the form 

T(r) = T Q (r) + (as)- 2 T 2 (r)-{-(as)- i Ur) + . . .. (6.2) 

After substituting (6.2) into (2.8), expanding cot (r/a) 
in a power series in r/a and equating the coefficients 
of like powers of (as)~ 2 , he obtained a system of 
simultaneous differential equations, which is written 
in terms of the operator 



as follows: 



L ^H& + li +1 ) T ' 



L(T )=0, 



Tfrp\ X cLIq , X dl 2 

Mi4J ~45 dx + 3 dx' 



2x 5 dT , x z dT 2 . x dT 4 



^^S^Tz+^V"^ 



3 dx 



(6.3a) 

(6.3b) 

(6.3c) 
(6.3d) 



where x=sr. 

Pekeris wrote the solutions of (6.3a) and (6.3b) 
in terms of the Hankel functions (corresponding to 
an implicit time function e i<at ) 



rr (sr) 2 



12 



m»(sr), 



and used 



(as)~ 2 T 2 (r) = -(lll2)(r/a) 2 H^(sr) 

as the term to estimate the error in using T as the 
solution of (2.8). 

We now proceed to extend Pekeris' treatment to 
solve for T 2n (r) for arbitrary n. The solution will be 
obtained in terms of cylinder functions of corre- 
spondingly higher order. Hence it is possible to 
obtain T(r) in the form (6.2) to any desired degree 
of accuracy for any distance r. 

Let Z n (x) be any cylinder function of order n. 
Define 

C n (x)=x n Z n (x) 
and 

L(T)=T"+x~ 1 T'+T, 

where primes indicate derivatives with respect to x. 
Then C n satisfies the differential equation. 



L(C n ) = (2nlx)Ci 



(6.4) 



Now a solution of (6.3a) is 

n' -7.-2/V 

^0 — ^ ^2) 



(6.5a) 



and since 
(6.3b) becomes 

L(T 2 )=-(\ft)x- l C 2 =-(l/\2^x- l Ci 

Hence from (6.4) it follows that 

T 2 = -(1/12) C 2 = - (1/12) (sr) 2 Z 2 (sr). (6.5b) 

Similar procedures may be applied to (6.3c), (6.3d), 
... by applying the following reduction formula to 
reduce the terms on the right-hand side to a sum 
of terms of the form 2nx~ ] CZ\ 



X C m — J-jC a ,k,m^m+k+a, k — 1,2, . . ., (6.6) 

a = 

where 

Oo,k,o =z 0? 

ali/c— a)\{m-\-a — 1)! 

(6.7) 

Equation (6.6) may be derived from the recurrence 
formulas for the cylinder functions. The solutions 
for T 4 and T 6 are given below: 



r4 ~~30 (78+ 160 C ' 4 ' 



?V— 



17 

Ml) 



C 4 - 



19 



Or 
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2520 b 120,960 



(6.5c) 



C,. (6.5d) 



In most cases of propagation around the earth, one 
or two terms of (6.2) give sufficient accuracy. 

Finally, the treatment may be extended in a simi- 
Lar way to the case of d^/d<£^0. Equation (2.4) 
then is replaced by the associated Legendre equation, 
for which we have obtained a solution in terms of 
generalized C functions 

^n, m\X) =Z X Zj m + n \X) . 

The details of this development will be reported 
elsewhere. 

7. Conclusions 

By a refinement of the procedure used in the usual 
earth-flattening approximation, we have shown that 
the problem of propagation around a spherical earth 
can be reduced to an exact equation of the same 
form. Thereby the earth-flattening procedure be- 
comes applicable to arbitrarily large heights and dis- 
tances. It is also found that existing solutions of the 
approximate equations can be re-evaluated to yield 
the exact solutions for slightly different refractive 
index distributions. 
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